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SUMMARY 


A  mathematioal  treatment  of  a  family  of  3id.table  designs  is  presented 
in  the  form  of  charts  shov/ing  the  inter-relation  of  density  of  loading,  sliver 
fraction  and  web  thickness. 
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I  nmoDucTioN 

Designing  solid  propellent  charges  is  an  art.  If  this  is  recognised 
then  the  methods  now  presented  will  be  accepted  in  their  right  perspective 
as  another  attempt  to  formulate  a  rational  basis  for  designing  charges  in 
which  burning  occurs  at  approximately  constant  pressure  on  the  surface  of  a 
central  conduit  of  star  form, 

The  condition  of  constant  burning  pressure  demands  constant  surface 
area  of  burning  except  for  the  initial  period  in  which  erosive  burning  may 
occur. 


Already  several  reports  have  been  written  on  the  mathematical  approach 
to  star-centre  charge  design.  Two  of  the  earliest  reports  known  to  the  writer 
are  those  of  references  1  and  2,  whilst  later  reports  giving  more  arithmetical 
results  are  listed  as  references  3,  4  and  5^  The  treatment  presented  in  this 
paper  is  new  and  the  results  are  in  a  convenient  form  for  assessment  work, 

2  EXPLANATION  OF  NOMENCLATURE 


A  typical  form  of  star-oentre  charge  is  shown  in  Fig.l  end  some  of  the 
defining  quantities  are  shewn  in  Pig, 2,  These,  and  others  not  capable  of 
pictorial  presentation,  are  given  below. 

The  mathematioal  exercise  is  feasible  only  if  the  stirfaoe  area  of 
burning  is  confined  to  the  conduit;  that  this  surface  is  the  product  of  the 
axial  length  of  the  conduit  and  the  perimeter  of  the  conduit  in  the  plane 
normal  to  the  axis  of  the  conduit,  and  that  the  axial  length  is  constant. 
Burning  at  constant  pressure,  now  defined  as  "neutral  burning",  therefore 
implies  that  the  perimeter  of  the  conduit  shall  be  constant. 

List  of  Symbols 

P’  Imaginary  initial  perimeter  (see  Pig, 2) 

P^  Actual  initial  perimeter 

P  Any  intermediate  perimeter 

P  Perimeter  at  instant  of  web  bum  out 
w 

R  Outermost  radius  of  charge  (inner  radius  of  motor  case  if  the  charge 
is  case-bonded) 

r^  Initial  radius  of  fillet  at  jumotion  of  lobe  and  web  (see  Pig,l) 

Tg  Initial  radius  of  lobe  tips 

N  Number  of  lobes 

w’  Web  thickness  corresponding  to  imaginary  initial  perimeter 

w  Web  thickness  corresponding  to  actual  initial  perimeter 

a  Constant,  where  P  =  2aR 

a*  Value  of  a  for  simultaneous  bum  out  of  web  and  lobe 

b'  The  distance  measured  normal  to  the  lobe  face  from  the  imaginary  initial 

jvinotion  of  the  lobe  face  and  the  web  to  the  oentre  line  of  the  lobe 
(see  Pig. 2) 
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’d  Specified  as  b'  -  Pj  (see  Fig. 2) 

P  Segment  angle  (including  semi-’web  and  semi-lobe  =  V^N 
^  Design  angle  (when  2^  includes  one  imaginary  web  portion) 

a  Design  angle  (see  Fig, 2) 


a*  Angle  subtended  by  circular'  arc  of  a  semi-sliver  at  the  imaginary 

initial  Junction  of  lobe  and  web,  at  instant  of  web  bum -out,  for 
the  case  of  prior  lobe  bum  through  (see  Fig,  10a) 


e 

X 

X' 


Sliver  fraction,  i, e. 


Total  cross-sectional  area  of  slivers  at  web  b\im~out 


•  Initial  cross-sectional  ai'ea  of  charge 

Density  of  loading,  i.e.  ^2iiial^o^s-s_ectiona^^ 

% 

Density  of  loading  for  the  imaginary  initial  perimeter 


3  THE  DESIGN  OF  STAR-CENTRE  CHARGES 
3. 1  Typical  charge 

The  cross-section  of  a  typical  star-centre  charge  is  shovm  in  Fig, 1, 
and  it  will  be  seen  that  there  are  no  sharp  comers.  The  apices  of  the  lobes 
are  rounded  to  redxice  the  initial  b\iming  perimeter  and  so  reduce  or  elimin¬ 
ate  the  high  initial  pressure  peaJc  that  would  otherwise  occur  (an  alternative 
treatment  is  suggested  in  reference  2),  Fillets  are  provided  at  the  junctions 
of  the  lobes  and  the  webs  to  eliminate  stress  concentrations,  which  might  lead 
to  cracking  of  the  charge  during  storage  or  firing. 


A  single  charge  segment  is  shov/n  in  Fig,  2,  together  with  generating 
lines  and  details  of  angles.  The  geometry  of  this  type  of  charge  design  is 
investigated  in  the  Appendix, 


3*2  The  effect  of  the  number  of  lobes 


The  final  perimeter  of  the  conduit  at  web  b\im-out  is  likely  to  be  only 
slightly  greater  than  2xR,  so  the  first  designs  investigated  were  those  giving 
neutral  burning  for  the  greater  part  of  the  effective  btirning  time,  and  having 
a  conduit  perimeter  of  27cR,  The  condition  for  neutral  bviming  is  derived  in 
the  Appendix  and  it  \’n.ll  be  seen  that  this  is  satisfied  by/ a  certain  value  of  the 
angle  (a  +  4>)>  this  hais  been  found  for  values  of  N  from  2  to  l6  and  the  results 
are  given  in  Table  I,  For  N  less  than  or  equal  to  Table  I  shovrs  that  the 
lobe  side  is  re-entrant  for  neutral  burning,  that  is  (a  +  ^)  is  greater  than 
90  degrees,  which  means  that, for  small  values  of  <f>j  and  sharp  comers,  an 
in^ossible  design  could  result  with  the  apices  overlapping.  The  maximum 
density  of  loading,  the  ratio  w/R,  and  the  sliver  fraction,  have  been  calcu¬ 
lated  for  the  above  condition  (P  =  27cR)  ,  and  are  tabulated  in  Table  II  for  N 
from  2  to  16,  for  designs  having  sharp  comers.  These  designs,  for  values  of 
N  of  5  or  less,  have  the  lobe  apices  just  meeting;  the  effect  of  rounding 
off  the  comers  is  dealt  with  in  the  Appendix, 

It  will  be  seen  that  the  maximum  density  of  loading  occurs  at  N  equal 
to  6,  with  only  a  slightly  lower  value  at  N  equal  to  5«  This  is  not  z’eally 
surprising,  since  consideration  of  the  angle  (oL  +  <j))  would  show  that  a 
maximum  density  of  loading  equal  to  1,0  yrould  be  achieved  when  (a  +  ^)  is 
equal  to  90  degrees;  for  neutral  bviming,  this  would  occur  with  a  value  of  N 
5*55  which  is  not  possible. 
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li.  DESIGN  CHANTS 
2f*1  General 

Detail  charts  have  "been  prodviced  for  star-centre  designs,  to  have 
neutral  huming  characteristics  (P  =  27cR),  for  values  of  N  of  5>  6,  7  and  8, 
and  these  are  shown  in  Pig.  3,  k>  5  and  6,  Basically,  the  density  of  loading 
is  plotted  against  weh  thickness/ charge  radius  (w/r)  and  the  radius  at  the 
-outer  star  point/charge  radius  (r^/R).  On  this  are  superin:5>osed  lines  of 
constant  sliver  area/motor  area  (Xe),  v/hich  occur  as  straight  lines,  and 
lines  of  constant  proportions  of  neutral  huming  (h/w).  Also  shown  in  a 
small  picture  is  the  ratio  of  the  final  perimeter  at  v/eh  hum-out  to  the 
normal  perimeter  (P^2'3tR)  plotted  against  the  design  angle  the  value  of 
i>  can  he  fouiid  from  the  main  chart,  since  the  values  are  shown  on  the  hound- 
aries  and  a  line  of  constant  ^  wovild  appear  as  a  strai^t  line  parallel  to 
a  line  of  constant  Xe» 

In  Fig.  3>  for  a  5~point  star  design, it  will  he  seen  that  the  upper 
part  of  the  diagram  has  heen  cut  off  at  a  density  of  loading  of  about  0,95« 

This  is  due  to  the  fvindamental  requirement  that  the  lohe  apices  cannot  overlap. 

The  effect  of  rounding  the  lohe  apex  to  reduce  the  initial  perimeter 
will  only  reduce  the  density  of  loading  X,  and  will  not  affect  the  values  of 
Xe,  or  P^/2xR.  The  proportion  of  neutral  huming  will  he  reduced  hy  r2/w 
at  the  beginning,  when  the  huming  is  progressive. 

Since  the  lines  of  constant  sliver/motor  area  (Xe)  are  not  evenly  spaced 
across  the  diagrams,  curves  of  sliver/motor  area  (Xe)  against  the  parameter 
(w  +  r^)/R  are  shovm  in  Pig.  7  for  values  of  N  of  5»  6,  7  and  8, 

For  given  weh  thickness/charge  radius  (w/R)  and  density  of  loading  X, 
the  sliver/motor  area  (Xe)  is  very  nearly  independent  of  the  nimiber  of  star 
points.  The  proportion  of  neutral  huming  (h/w)  reduces,  and  the  final 
perimeter  increases,  with  increase  in  the  number  of  star  points.  Also,  for 
a  given  weh  thickness/charge  radius  (w/r),  the  range  of  density  of  loading 
tha.t  can  he  considered  is  reduced  as  N  is  increased,  the  maxitnum  heiug 
reduced  and  the  minimum  increased.  Prom  these  points  of  view,  it  is  suggested 
that  5  or  6  pointed  star-centre  designs  are  likely  to  he  the  most  usefvil. 

1|..2  Six-pointed  star-centre  designs  v/ith  values  of  a  other  than  % 

It  was  considered  that  some  investigation  should  he  made  into  the  effect 
of  veirying  the  constant  a.  In  Pig.  8,  the  density  of  loading  with  sharp 
comers  (X*)  is  shovm  as  a  fimction  of  weh  thickness/charge  radius  and  the 
constant  a.  Also  shown  are  lines  of  constant  sliver  area/motor  area  and 
lines  of  constant  final  perimeter/normal  perimeter  (P■^Y/2aR).  The  value  of 
the  design  angle  <f>  is  shown  in  Pig.  9  against  v/eh  thickness/chai'ge  radius 
and  the  constant  a.  With  sharp  comers,  lines  of  constant  proportions  of 
neutral  huming  (h ’/■«■’)  are  also  shown. 

As  is  shown  in  the  Appendix,  for  a  greater  than  a*,  the  vreh  will  hum 
out  before  the  star  side  and  there  will  he  neutral  huming  until  weh  hum-out. 
In  the  present  case  a*  is  slightly  less  than 

Let  \as  now  consider  the  effects  of  varying  the  value  of  a.  Suppose 
we  start  at  a  lairge  value  of  a ,  At  constant  density  of  loading,  as  a  is 
reduced  we  have  the  following  effects; 
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(i)  the  sliver  fraction  is  reduced,  since  X  is  constant, 


(ii)  web  thickness/charge  radius  is  increased, 


lere  is  no  change  in  the  final  perimeter/norraal  perimeter 
until  a  =  a*,  but  as  a  is  further  reduced,  this  ratio 

increases, 

IS  an  alternative  to  (iii)  the  proportion  of  neutral  burning 
reduces  as  a  is  reduced  below  a*. 


Due  to  the  method  of  filling  rocket  motor  tubes  with  plastic  propellent 
some  increase  in  perimeter  can  be  permitted  if  the  burning  surface  is  to  be 
kept  constant,  since  the  initial  length  of  the  charge  is  greater  than  the 
final  length,  as  the  propellent  is  pressed  into  the  head  and  rear  danes  as 
well  as  the  parallel  portion  of  the  tube.  Since  the  reduction  in  length  is 
\mlikely  to  exceed  10^,  unless  the  charge  length/diameter  is  low,  the  maximum 
value  of  final  perimeter/normal  perimeter  vrould  be  1.10.  Also,  there  do  not 
appear  to  be  any  very  valid  reasons  for  having  a  greater  than  a"^,  so 
that  the  reasonable  range  of  a  is  from  2.9  to 

4«3  Approximate  equations  for  six-pointed  sta3>-centre  designs 

For  project  assessment  v/ork,  it  was  considered  that  it  would  be  useful 
to  have  sane  approximate  equations  relating  density  of  loading,  web  thickness, 
sliver,  etc.  From  the  earlier  discussion,  several  reasons  have  been  given 
for  using  5  or  6-pointed  star  centres,  and  for  this  case  the  6-pointed  star 
was  used,  with  a  perimeter  of  27dl, 

Assuming  that  a  density  of  loading  less  than  0. 6  would  not  be  of 
interest,  the  following  equations  were  obtained,  and  the  R.M.S.  error  is 
also  given; 

As  linear  equations  R,M,S,  error 


X‘  =  2.661888  (w»/r)  -  0.167889  , 
Xs  =  0.658981  (w’/R)  -  0.168513  , 


0. 001^722 


0.00[f791 


X’-Xe  =  2.002907  (w‘/R)  +  0.000622^  , 


0.000069 


w*/R  =  0.009055  ^  +  0.428325  for  ^  in  degrees. 


0.003530 


R.M.S.  error 


As  parabolic  equations 


X'  =  3.363424  (w’/R)^  +  0.254098  (w«/R)  +  0.257747  , 

Xe  =  3.412416  (w«/R)^  -  1.783882  (w'/R)  +  0.263323  , 

X’-Xe  =  -  0.048993  (w’/R)^  +  2.037979  (w’/R)  -  0.005576  , 


0.000752 


0.000766 


0.000014 


w’/R  =  -  0.000211  0,005999  ^+0.421195  for  0  in  degrees. 


0,000191 


Now,  as  shown  in  the  Appendix,  the  actual  density  of  loading  is  given 

by: 
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X  =  X'  -  (2r/R)  -  ir^Rf 
and  (w'/K)  =  (w/R)  +  (r^/R) 

and  =  (1  -  P^/P)^ 


e.g.  it  may  he  found  that  for  a  given  propellent  under  certain  conditions,  a 
value  of  P./p  of  0,9  will  he  stiff icient  to  remove  an  initial  pressure  peak  and 
a  value  of  0,005  for  r^/R  v/ill  he  adequate  to  safeguard  against  charge  cracking. 
Then  - ' 


X  =  X'  -  0.02 

(w/R)  =  (w'/R)  -  0.005 

and 

X  -  Xe  =  (X*  -  Xe)  -  0.02  . 


It  will  he  realised  that  (X  -  Xe)  represents  the  charge  area  humt  away 
at  hum-out,  and  the  equations  given  represent  extremely  good  approximations 
to  this. 
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TJffiLE  I 

Conditions  fctr  neutral  'bumlnK 


N 

1 

P 

u  +  ^ 

Radian 

Degree,  minute 

Radian 

Degree,  minute 

2 

1.57079633 

90° 

2.79838605 

160°  20.143’ 

3 

1.0^719755 

60° 

2.18957409 

125°  27.201’ 

4 

0.78539816 

45° 

1.86372955 

106°  47.030’ 

5 

0, 62831853 

36° 

1.65580688 

94°  52.245’ 

6 

0.52359878 

30° 

1.50921626 

86°  28,303* 

7 

0.44879895 

25°  42.857’ 

1.39901178 

80°  9.448’ 

8 

O.392699O8 

•• 

0 

0 

CM 

CM 

1,31237148 

75°  11.601’ 

9 

0.34906585 

20° 

1,24197868 

71°  9.608’ 

10 

0.31415927 

18° 

1.18332827  I 

67°  47.983’ 

12 

0.26179939 

15° 

1.09042934  1 

62°  28.620’ 

14 

0.22439948 

12°  51.429' 

1.01940817  j 

58°  24.467’ 

16 

1  0.19634954 

1 - -  - 

11°  15’ 

i 

0.96279471  j 

55°  9.844’  i 

TABLE  II 

Maximum  density  of  loading  for  neutral  bimilnp:  and  P  =  2?^ 


N 

w’/R 

X’ 

2 

0.294845 

0.611400 

3 

0.353094 

0.752720 

4 

0.392545 

0.862372 

5 

0.420910 

0.954628 

6 

0.2^21582 

0.964403 

7 

I 

0.398500 

0.897179 

8 

0.378057 

0.841057 

9  ! 

j  0.359935 

0.793307 

10 

0.343798 

0.752057 

12 

0.316331 

0.684073 

14 

0.293810 

0.630a^8 

16 

!  0.274970  ! 

0.585842 

1 


Xe 

Pw'/2xR 

0.020985 

1,021495 

0.047644 

1.022683 

0.075791 

1.029316 

0.111233 

1.029335 

0.119542 

1.030276 

0.098275 

1.032533 

0.082943 

1.033645 

0.071407 

1.034093 

0.062442 

1,034144 

0.049475 

1.033608 

0.aH0608 

1.032685 

0.034205 

1.031634 
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2  lu 

Some  consideration  has  already  "been  given  to  this  type  of  design 

but  with  various  limitations,  such  as  ^  =  0,  r2  =  0,  or  r^/R  =  0,06,  To 

simplify  the  equations  we  will  assume  initially  that  r-  =  0  and  use  w’ 
instead  of  (w  +  r^). 

The  initial  perimeter  P’  is  given  by 


£1 

2N 


(R  -  w’)  ^  + 


(R-wQ  sin  (P-  0) 

COS  (a  +  ^  -  (3) 


(1) 


and  after  a  distance  x  has  burnt  the  perimeter  P  is  given  by 


P 

X 

2N 


(R-w> 


/ 


y 


+x)  (f)  +  xa  + 


y  . 

(R  -w*  )  sin  (P  -  0) 

cos  (a  +  ^  -  3) 


-  X  tan  (ct  +  ^  -  3) 


(2) 


provided  that  the  star  side  has  not  bvimt  out. 

Until  the  star  side  or  the  v/eb  bvims  out,  whichever  occurs  first,  v/e 
have,  by  differentiating  (2) 


JL 

2N 


a,  +  <f>  -  tan  (a  +  ^  -  3) 


(3) 


Density  of  loading 

The  cross-sectional  area  of  the  charge  (x  R  X*)  can  be  shown  to  be 

N  S  _  (R_v,')^  #  -  .(P.  -  ,(P  rt,.C93_a' 

COS  (a  +  ^  —  3) 


or 


X> 


-i  A  ^  -  /-i  sin  (3~s^)  003  a  /,  n 

“  ^  “  V  V  V  y  ^  cos  (a4.^i-3)  • 


Sliver  area 


The  portion  of  the  charge  referred  to  as  sliver  is  shcrwn  in  Pig,  10a 
and  Pig,  10b,  where  it  v/ill  be  seen  that  two  equations  are  required,  according 
to  v/hether  the  star  side  b\ims  out  before  the  vreb  or  the  converse.  If  the 
side  has  burnt  out  first  as  in  Pig,  10a 
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Sliver  area/N 


R^(P  -  0)  -  (w*)^  a’  -  w’(R  -  w’)  sin  a’ 


where 


a* 


(P  -  0)  +  sin 


[(R  -  w’)  sin  (P  -  0)  w’]  j 


or  sliver  area/motor  area, 


Xe 


=  1 


/w’\^  a’  /w’\/,  w’\  sin  a’ 

f-{j)  T  - -  t) -J- • 


(5) 


For  the  case  shovm  in  Fig, lOb,  where  the  web  bums  out  before  the  star- 
side,  it  can  be  shown  that 


Sliver  area/N  =  R^(P  -  <f>)  +  (w’)^  [tan  (a  +  ^  -  P)  -  a] 

_  (R  -  w')  sin  (P  -  a  ((g  _  ^  2  w'i  , 

COS  (oc  +  ^  -  P) 


or  sliver  area/motor  area. 


^  ^  p  +  i  r]  P 


-  P)  -> 


-[( 


P  cos  (a  +  ^  -  P) 


sin  (P  -  ^)  r/ 


w' 


2w’ 


1  -  Tj  cos  «  +  R 


J 


.  (6) 


Effect  of  rounding  outer  star  point  i,e,  r^  =f=  0 


Decrease  in  density  of  loading  =  -  AX^ 

r . 


•1 


TtR 


P  dx 
X 


- 1  -  (t) 


-  11  - 

CONFIDENTIAL 


CONFIDENTIAL 


Technical  Note  No»  190 
Appendix 


Effect  of  rounding  lobe  apex  i*e.  Tg  +  0 


Reduction  in  initial  perimeter 


=  2N  r2  [tan  (a  +  ^  -  (3)  -  (a  +  ,  (8) 

Reduction  in  density  of  loading  =  -  AXg 


[tan  (a  +  ^  -  3)“(a  +  ^  -  P)]  • 


(9) 


Final  perimeter 

If  the  star  side  bums  out  first,  the  final  perimeter  at  v/eb  biim-out, 
is  given  by 

Py2N  =  R(P  -  <f>)  +  w’a'  .  (10) 


If  the  star  side  has  not  bximt  out 


P 


w 

2N 


R^  +  vr'a  + 


(R  w*)  sin  (P  - 

cos  (a  +  ^  -  P) 


w‘  tan  (ot  +  ^  -  P) 


(11) 


Neutral  burning 


In  this  case,  it  is  necessary  to  make  dP^dx  =  0,  or 


(a  +  0)  =  tan  (a  +  ^  -  P) 


(12) 


where 


P  =  Vn  • 

Equation  (l2)  has  been  solved  by  iteration  for  N  =  2,  3>4»«*9>  10,  12,  Ilf,  1 
and  the  values  of  (a  +  0)  are  given  in  Table  I  in  both  radian  and  circular 
measure# 

For  a  given  value  of  N,  increasing  the  value  of  (a  +  ^)  from  that  for 
neutral  burning  will  make  dP^dx  negative,  and  reducing  the  value  of  (a  +  0) 

will  give  progressive  bximing. 

Initial  perimeter  =  ?aR 

Substituting  this  value  in  equation  (l)  we  obtain 


-  12  - 

CONFIDENTIAL 


CONFIDENTIAL 


Technical  Note  No»  190 
Appendix 


+  sitx  (P  >>  0)  _  ~j  ^ 
cos  (a  +  0  -  3)J 


(13) 


Now  (a  +  f)  and  3  are  constants  for  given  N,  and  equation  (l3)  enables  the 
value  of  (l  -  w’/r)  to  be  found  for  various  values  of  0. 


Effect  of  rovmding  outer  star  points  end  lobe  apices 

In  this  cane,  when  r^  is  small,  is  constant  at  2aR,  Hence 
equations  (7)  and  (9)  become 

-  =  2a  r^/%R  (7a) 

-  AX2  = 


Density  of  loading.  X 

This  quantity  is  foxmd  from  equation  (4)  using  the  value  of  w‘/R  from 
equation  (13),  and  equations  (7a)  and  (9a). 

Sliver  fraction 


The  value  of  sliver  fraction  is  obtained  by  using  equation  (5)  or  (6) 
together  with  the  density  of  loading  found  above.  The  value  of  a  will 
determine  whether  equation  (5)  or  equation  (6)  should  be  used.  The  star  side 
will  bum  out  at  the  same  time  as  the  web,  if  the  length  b'  in  Pig. 2  is  equal 
to  w’.  That  is,  if 

(R  -  w’)  sin  (3-0)  =  w'  sin  (a  +  0  -  3)  > 


and  this  will  be  satisfied  -vrfien 


a  =  a* 


N(a  +  0) 


sin  (3  -  0)  +  0  cos  (g  +0-3) 

sin  (3  -  0)  +  sin  (a  +  0  -  3) 


(14) 


The  value  of  a*  has  been  evaluated  for  N  between  2  and  I6  and  0 
between  0  and  3  and  the  results  are  shown  in  Pig, 11.  Also  shown  in  this 
figure  is  a  line  for  a  =  90  deg.  which  occxxrs  when  the  star  points  just 
meet  for  r^  and  =  zero.  If  a  is  less  than  a* ,  the  side  vnll  bum  out 

before  the  web,  and  if  a  is  greater  than  a*  ,  the  web  will  bum  out  first. 
Proportion  of  neutral  burning 

This  is  the  ratio  (b/w)  or  (b’  -  r^)/(w’  -  r^). 
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Now 


^  LR  ~  wQ  sinJP  -  j6l 
sin  (a  +  0  -  13)  * 


and,  therefore, 


h 

w 


^(R  -  w')  sin  (P  -  0)  -  r^  sin  («  +  0  -  (3) 
(w*  -  r^)  sin  («  +  0  -  (3) 


(15) 
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FIG.  5  NEUTRAL  BURNING  STAR-CENTRE  CHARGE  DESIGN 
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